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Abstract 

Spintronics is the science and technology of electric control over spin currents 
in solid-state-based devices. Recent advances have demonstrated a coupling 
between electronic spin currents in non-magnetic metals and magnons in 
magnetic insulators. The coupling is due to spin transfer and spin pumping 
at interfaces between the normal metals and magnetic insulators. In this 
Chapter, we review these developments and the prospects they raise for 
electric control of quasi-equilibrium magnon Bose-Einstein condensates and 
spin superfluidity. 


1.1 Introduction 

In undergraduate texts, Bose-Einstein condensation is a phase transition 
that occurs in an ideal Bose gas at large enough densities. In developing set¬ 
ups to observe Bose-Einstein condensation in a way analogous to this text¬ 
book treatment, one faces several challenges. First of all, conserved bosons 
in condensed-matter systems are composite particles. For example, Cooper 
pairs of electrons are bosons but the Bardeen-Cooper-Schieffer regime of 
condensation (in which the electrons condense as a result of a weak and 
attractive effective interaction) is rather different from the physics of non¬ 
interacting point-like bosons [T]. The size of the Cooper pairs, determined in 
the weakly-interacting limit by the coherence length ^ ~ ep/kpA with kp 
the Fermi wave number, ep the Fermi energy, and A the superconducting 
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gap, is much larger than the distance between the pairs. This latter dis¬ 
tance is estimated from the fraction A/ep of the electrons that form Cooper 
pairs, so that the pair density is Up ~ kpA/ep ~ kpj^, and we have that 
~ > 1 - 

Cold atoms [2], on the other hand, are composite bosons where the in¬ 
ternal degrees of freedom are typically at much higher energies than their 
temperature, so that they can be considered as point particles. The bosonic 
atoms in magnetically-trapped ultracold atomic vapors are, however, not 
conserved and have a finite lifetime because atoms may escape from the 
trap (both as single atoms or after molecule-forming collision processes). 
The same mechanism also prevents the system from reaching its true ther¬ 
modynamic ground state (a state that is most likely a solid). Relaxing the 
requirement of strict boson conservation, one may consider Bose-Einstein 
condensation as a quasi-equilibrium phenomenon. Prime examples of quasi¬ 
equilibrium Bose-Einstein condensation of non-conserved particles are con¬ 
densates of photons [3] and (exciton-)polaritons [UIl]. In true equilibrium, 
these particles would not be described by a Bose-Einstein distribution func¬ 
tion with nonzero chemical potential and for both photons and polaritons 
the system is maintained in quasi-equilibrinm at nonzero chemical potential 
by external pumping. The mechanism of reaching quasi-equilibrium is very 
different for each of these systems. 

In this Chapter we focus on magnons, quasi-particles corresponding to the 
quantized fluctuations of the magnetic order parameter, in magnetic insula¬ 
tors. These can undergo Bose-Einstein condensation as an equilibrium PEj 
or a quasi-equilibrium [8] phenomenon. Here, we focus on quasi-equilibrium 
condensation of magnons. This is the result of excited magnons coming to 
quasi-equilibrium by magnon-conserving scattering processes, before they 
dissipate energy to the lattice and relax. In experiments BM on quasi¬ 
equilibrium magnon condensation in solid-state magnetic insulators the ex¬ 
citation of magnons is achieved by microwave pumping. Recent advances in 
spintronics, however, have demonstrated interactions between magnons and 
electrons at interfaces between metals and magnetic insulators |in) . As we 
will discuss in detail in this Chapter, this opens the possibility of direct- 
current (DC) pumping of the magnetic system (see Fig. |1.1[ ), and to achieve 
quasi-eqnilibrium magnon Bose-Einstein condensation in a solid-state DC 
transport experiment. On top of this, the integration of quasi-equilibrium 
magnon Bose-Einstein condensation with electronics opens the possibility 
for electronic transport probes of the associated spin snperfluidity. 

In the remainder of this Chapter we first discuss the recent developments 
in spintronics that open the possibility of manipulating the magnetization in 
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Figure 1.1 Schematic of a normal-metal/magnetic-insulator heterostruc¬ 
ture, in which a finite chemical potential of magnons can be induced by a 
nonequilibrium electron-spin accumulation fi. In the idealized case with no 
magnon-number relaxation in the insulator, the metallic reservoir acts as 
a bath supplying magnons and energy into the ferromagnet, with temper¬ 
ature T and chemical potential p. SSn is the entropic change in the metal 
bath associated with the transfer of energy 6 U and creation of SN magnons 
in the ferromagnet (relative to the ordered spin orientation n). 



magnetic insulators by means of electrical current. Hereafter, we discuss the 
phase diagram of the DC-pumped quasi-equilibrium magnon gas. Finally, we 
briefly discuss signatures of the resulting spin superfluidity and give future 
prospects. 


1.2 Spintronics 

In this section we outline recent discoveries in spintronics that enable the 
integration of quasi-equilibrium magnon Bose-Einstein condensation with 
the control of spin currents. The first of these phenomena is the existence 
of spin currents across an interface between a normal metal and a magnetic 
insulator. The second concerns the generation of spin currents from charge 
current (or vice versa) via the (inverse) spin Hall effect. 


1.2.1 Interfacial spin currents 

We consider the set-up in Fig. |1.2| a magnetic insulator sandwiched between 
two normal metals. This geometry is chosen for convenience as it makes the 
theoretical analysis of the spin transport effectively one-dimensional. In an 
experiment it is easiest to put the metals on top of a thin film of a magnetic 
insulator. In terms of materials, an often studied system is the magnetic in¬ 
sulator Ytrrium-Iron-Garnett (YIG) interfaced with the non-magnetic nor- 
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Figure 1.2 A set-up that we consider in this Chapter. A quasi-equilibrated 
magnetic insulator at temperature and chemical potential fj,m is con¬ 
nected to left and right by metallic reservoirs with spin accumulations fiL 
and flu, and temperatures and Tr, respectively. The spin accumulation 
in the reservoirs results from the spin Hall effect in response to currents in 
the j/-direction (shown for left reservoir) and/or injected spin current from 
the magnetic insulator. Any spin current injected from the insulator to the 
reservoirs results in an inverse spin Hall voltage which can be measured 
(shown for right reservoir). The spin quantization axis is defined by the 
external field that we choose in the z-direction so that the equilibrium spin 
density points in the —z-direction. 


mal metal platinum. We assume that the electrons in the normal metal 
on the left (right) have an electron spin accumulation that is nonzero at 
the normal-metal side of the interface between normal metal and magnetic 
insulator and is denoted by This spin accumulation is an out- 

of-equilibrium electrochemical potential imbalance between up and down 
electrons, fiL = fj.^ — m (with a similar dehnition of fj-ji). While in a more 
general set-up, the spin accumulation would be a vector, here we assume 
that only its z-component is hnite, as illustrated by the electron spins in 
Fig. 1.2 in the normal metal. How the spin accumulation is established is 


discussed below. We further assume that magnons in the magnetic insulator 
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are in quasi-equilibrium characterized by a magnon chemical potential pLm 
and temperature Tm- The metallic leads are thermally biased, with Tp and 
Tpi for the left and right reservoirs, respectively. 

We first focus on the left interface. At the interface, there exists an inter¬ 
face exchange coupling between the localized spins in the magnetic insulator, 
and the itinerant spins in the metal. Assuming this coupling is isotropic, a 
phenomenological expression is 

Vint = J dxdx'F(x,x')S(x) • s(x') , 


where S(x) corresponds to the spin density of the localized spins in the 
insulator, and s(x) the spin density of the electrons in the metal. The latter 
is given by 


s(x) 


I 

cr,(T'e{td} 


in terms of the Pauli matrices r and electron creation and annihilation oper¬ 
ators V’ct(x) and 'i/’o-(x). Finally, the matrix elements P(x,x') decay rapidly 
away from the interface. 

We first assume that the spin density in the magnetic insulator can be 
treated classically and is homogeneous, and, moreover, that the temperature 
is relatively low so that we can approximate (S) ~ hsn, with the density 
s = S/v where S is the total spin per unit cell, v the volume of a unit cell in 
the magnetic insulator, and n a unit vector in the direction of spin density. 
The ^-component of the spin current from the magnetic insulator into the 
left reservoir, per unit area, is then given by m (we define positive spin 
current as flowing to the right) 


.nt _ ^9n 


( 1 . 1 ) 


with the so-called spin-mixing conductance (in units of m“^ and disre¬ 
garding its imaginary component). As we shall see, this latter conductance 
characterizes the efficiency of spin transport across the interface. It can be 
straightforwardly calculated, e.g., by using perturbation theory in the cou¬ 
pling l/(x, x'). A microscopic expression is not needed at this point. The 
mixing conductance can also be determined from ab initio calculations [1^ 
or from experiments [10]. For interfaces between YIG and Pt the mixing 
conductance is estimated to be up to 5 nm“^ depending on interface quality 
nans]. We note that the above expression captures the spin-pumping con¬ 
tribution and not the spin-transfer contribution that results from the spin 
accumulation in the normal metal. 
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We now consider the above expression for magnons. Using a linearized 
Holstein-Primakoff transformation, we have for the (circular) magnon anni¬ 
hilation operator that 


^ = \/2 - iSfiy) 


where we assumed, as in Fig. 1.2, the magnetic order to be in the —z- 


direction, so that n ~ {dhx, Shy, —1). Inserting this in the expression for the 
spin current leads to the replacement 


n X 


dn 

dt 


4 1 ^ 


with hcvk the magnon dispersion, and where we inserted an additional fac¬ 
tor of two to incorporate constructive interference of magnon modes at the 
interface [H]. Here, we have made a Fourier transform so that the number 
of magnons at momentum k is Uk = (^k^k), and V is the volume of the 
magnetic insulator. Furthermore, we normal-ordered the magnon creation 
and annihilation operators and have kept only the expectation values of ¥b 


as we are interested in thermal magnons. Inserting this result in Eq. (1.1) 
yields in first instance 


■int 
J s 


TTS V 


We now consider the situation of quasi-equilbrium magnons at chemical 
potential /Jm and temperature so that 

/huJk - hm\ 

keTm J ’ 

with nsix) = [e^ — 1]“^ the Bose-Einstein distribution function. In equilib¬ 
rium there is no spin current accross the interface. To account for this we 
generalize our treatment with the replacement (see Ref. |15j l 


riB 


hujk - k-m \ 
ksTm ) 


ns 


- hm \ 

ksTm ) 


- riB 


huJk - hL \ 

kBTi J ■ 


The first contribution corresponds to spin pumping, the emission of spin 
current from an excited magnet. The second term reflects spin transfer, 
the injection and absorption of spin current into a magnet. Over the last 
decade, both these phenomena have been investigated extensively for a clas¬ 
sical magnetization, i.e., not including contributions of magnons (see e.g. 
Ref. [16] and Ref. |lZ|)- Putting the above together, and rewriting momen¬ 
tum summations as energy integrations involving the magnon density of 
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states D(e), we finally find the result 

ir = -— fdeDie) {e - Lf 
ns J L V 

In linear response this yields 


ksTm 


- riB 


e- BL 

ksTp 


■int 

Js 


a 


int 


hA 


T int 

{BL - Bm) + - Tm) 


( 1 . 2 ) 


(1.3) 


which defines an interface spin conductivity cr®"'* = 3C(3/2)h5'|-^/27rsA^ and 
the coefficient L'^sse ~ ^^C{^/‘^)kB9^i/^nsA^. Here, we assumed a gapless 
quadratic magnon dispersion TiWk = Js^ in terms of the spin stiffness J*, so 
that the magnon density of states is D{e) = . Furthermore, we 

have introduced the length scale A = ^AnJgjkBTmi the thermal DeBroglie 
wave length for the magnons. Eqs. (1.2) and ([Tl are the main expressions 
for the magnon spin current accross the interface that we will use below. The 
interface between the magnetic insulator and normal metal on the right can 
be treated analogously. The contribution proportional to the temperature 
difference in Eq. (1.3) and determined by the coefficient L'^sse gives rise to an 


interface contribution to the so-called spin Seebeck effect, i.e., a spin current 
as a result of a temperature gradient [laiin]. This effect, and its reciprocal 
dubbed the spin Peltier effect [20], are subject of intense investigation |2T|. 


1.2.2 spin Hall effect 


Having addressed the spin current through the interface, we now turn to the 
question of what establishes the spin accumulations p l and pR in the normal 
metals on the left and right. These are the result of spin-current injection 
from the magnetic insulator into the normal metals, combined with spin 
current due to the spin Hall effect. The spin Hall effect arises due to spin- 
orbit coupling and is reflected as a spin current with a spin polarization and 
a spatial direction transverse to an applied electric field in a metal 


In the geometry of Eig. 1.2 the electric field (or charge current) is flowing in 
the y-direction, and gives rise to a spin current that flows in the x-direction 
with spin polarization in the 2 -direction. We now focus again on the left 
metallic reservoir and in first instance assume transport is diffusive with 
weak spin-orbit interactions (so that the scattering mean free path is shorter 
than the spin-flip diffusion length). Even though the diffusive approach is 
not generally applicable to experimental situations, as detailed below, it for 
now serves our pedagogical purpose. The set of equations that describes the 
coupled spin and charge dynamics in the normal metal are then given by 












(within our quasi-one-dimensional geometry) |24j 


jc = aE+ 5^/iL , 

2 e 

‘^js = - ctshE . (1.4) 

In the above the electric field E and charge current jc are in the y-direction 
and the electron charge is e. The charge conductivity a and spin Hall con¬ 
ductivity asH are both in units of m“^. The second term in the first 
equation is the charge current that results from a gradient in the spin accu¬ 
mulation via the inverse spin Hall effect [2S] . This is the Onsager reciprocal 
of the spin Hall effect and it is thus governed by the same coefficient ctsh- 
The inverse spin Hall effect is a powerful means to detect spin current elec¬ 
trically, as discussed in more detail below. Writing asH = Gsh(^, the spin 
Hall effect is quantified in terms of the dimensionless quantity (dubbed spin 
Hall angle) Osh- For Pt, Osh ~ 0.05 and for Ta its magnitude is similar, but 
the sign is opposite [26] . 

The above equations have to be complemented with a continuity-like equa¬ 
tion. 


djs _ ^ 

dx ’ 

where the rate per unit volume T phenomenologically expresses spin-flip 
relaxation in the metal on the left of the insulator. Insertion of the expression 
for the spin current into the latter equation yields the spin-diffusion equation 


dx^ P' ’ 


(1.5) 


with the spin-flip diffusion length of the left lead i = \fohJ\PE. 

We now assume the interface between the left metallic reservoir and the 
magnetic insulator is at position x = 0 and that the thickness of the metal 
is Lx in the x-direction and in the z-direction. The equations for the spin 
current and the spin accumulation are solved with the boundary conditions 
js(x = —Lx) = 0 and js(x = 0) = jf^^. Using the linear-response expression 
in Eq. ([Tl we find that the spin accumulation at the interface in the left 
reservoir is given by (assuming Osh ^ 1) 


Osnlch^E 

1 — cosh 

-h 2 eLL 2 [o-fVm + f^PsSEC^rn “ Tg)] cosh 

2eLxL,ai^t 

cosh(^) + A^sinh(^) 



where Ic = jcLxL^ is the total current through the normal metal (where the 
current density is assumed to be homogeneous). Note that the above result 
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Figure 1.3 Spin accumulations p,L and Pr, and magnon chemical poten¬ 
tial firm that build up in the metallic reservoirs and magnetic insulator 
in response to a thermal gradient across the magnetic insulator. The spin 
accumulations and magnon chemical potential are nonzero over distances 
~ i and ~ irn (in metals and insulator, respectively) away from the inter¬ 
faces. The jump across the interfaces is determined by the interface spin 
conductivity cr*"*. 


shows that in the limit of no spin relaxation, i.e., —)• oo, no net spin 

current is flowing across the interface as js is constant and zero (according 
to the boundary condition at —Lx). Again we mention for completeness that 
an analogous derivation holds for the right metallic reservoir. 


1.2.3 Inverse spin Hall voltage and the spin Seebeck effect 

Spin current can be detected electrically via the inverse spin Hall effect as 
spin current injected into the metallic reservoirs gives a inverse spin Hall 
voltage in the y-direction across the reservoirs (in an open circuit geom¬ 
etry). As a further illustration, we extend the drift-diffusion treatment to 
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capture the spin Seebeck effect. That is, we consider the inverse spin Hall 
voltage generated across the reservoirs as a result of a temperature gradient 
across the magnetic insulator. We assume, for simplicity, the temperature 
drops linearly Tm(x) = (Tr - TL)x/Lm + Tr, with Tr > Tr (see Fig. [Q] ). 
We do not consider a temperature gradient across the metallic reservoirs, 
nor do we consider the effect of the temperature difference across the inter¬ 
face (governed by the Kapitza resistance). In the diffusive limit of magnon 
transport, the magnetic insulator is described by similar equations as the 
diffusive normal metal, i.e., 


Js — 

dx"^ 


(Ts OfJjYn Y 

- ^SSE —^— 

OX 


% dx 
_ l^m 

~ w 


( 1 . 6 ) 


where (Ts is the spin conductivity of the magnetic insulator, Lsse its bulk 
spin Seebeck coefficient, and ^rn the thermal-magnon propagation length. 
Contrary to the case of Pt normal metals, a diffusive approach is expected to 
be appropriate for thermal magnons, albeit that the various microscopic de¬ 
tails (such as the role of complicated spin-wave dispersions, strong magnon- 
phonon and magnon-magnon interactions) are not well-understood and may 
lead to a non-trivial dependence of transport coefficients on magnetic field 
and temperature. The transport coefficients Lrse and Us and the length 
scale are therefore neither experimentally nor theoretically very well un¬ 
derstood for YIG at present. This, together with the complication of deter¬ 
mining magnon temperature profiles from experimentally applied tempera¬ 
ture differences |27ll28j . hinders a full and quantitative understanding of the 
spin Seebeck effect in YIG. Here our purpose is to discuss a simple example 
of diffusive spin transport that will be contrasted with the case of superfluid 
magnons later on. 

We solve the spin diffusion equations in the magnetic insulator subject 
to the boundary conditions that the spin current vanishes at the left and 
right boundary of the system. For both interfaces we use the boundary 


condition in Eq. (1.3) and we take zero charge current in both reservoirs. 


The properties of the metallic reservoirs, as well as their interfaces with 
the magnetic insulator, are taken equal. Fig. |1.3| shows a schematic plot of 
the spin accumulation and magnon chemical potential that build up in the 
metallic reservoirs and magnetic insulator as a result of the thermal gradient. 
Within a distance ~ I'm from the interfaces the magnon chemical potential is 
nonzero. Similarly, the spin accumulation in the metallic reservoirs is nonzero 
within a distance ~ The jump from spin accumulation to magnon chemical 
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potential across the interface is inversely proportional to the interface spin 
conductivity cr®"’*. 

The spin Seebeck coefficient is defined as SssE = Vishe/{Tl — Tr), with 
the inverse spin Hall voltage Vjshe found by computing the electric field via 


Eqs. (1.4) from the injected spin current from magnetic insulator to metallic 
reservoir, and averaging over the x-direction. We consider now the limit 
Lm which allows us to focus on one interface, which we choose to be 

the right one. Within the drift-diffusion theory the spin Seebeck coefficient 
is ultimately found as 


SsSE = — 


OsH^LyimLsSE 


‘ZcLxLyyi 


a 


int 


+ 


T)' 


-U — 


(1.7) 


where Ly is the length of the normal metals in the y-direction and where 
we also took the limit Lx ^ i- Within our drift-diffusion theory the spin 
Seebeck coefficient (times Lm) saturates as a function of Lm/^m to the value 
determined by the above result. This saturation results from the fact that 
only magnons within a length (.m from the interface contribute to the spin 
Seebeck voltage [29l |30] . 

In the above treatment the properties of the insulator are characterized by 
the phenomenological constants im, Lsse and as- The discussion is made 
more quantitative by computing these in the relaxation-time approxima¬ 
tion in which as rsj Jst/A^ and Lsse rsj JsksT/fih? , where r is the magnon 
transport mean free time. Furthermore, the magnon spin propagation length 
is then given by im ~ Vm\/TTsr, where Vm = 2\/ JgksT/h is the magnon 
thermal velocity and Tsr is the magnon spin-relaxation time. Here, r is the 
result of various magnon conserving and non-conserving relaxation mecha¬ 
nisms such as magnon-phonon scattering, magnon-magnon (Umklapp) scat¬ 
tering, and scattering of magnons with impurities. The spin-relaxation time 
Tsr acquires contributions only from processes that do not conserve magnon 
number. As we mentioned before, such relaxation and scattering processes 
are at present not fully understood. We remark that, by assuming Gilbert 
damping as the only relxation mechanism (which would be appropriate in 
clean systems at low temperatures), however, both time scales are on the 
order of h/aksT, with a the Gilbert damping constant (which for YIG can 
be as small as 10“^). 

Because the spin diffusion length of Pt (which is on the order of a few 
nm at room temperature [26]) is not large compared to the mean free path, 
and because in experiments the Pt layers may be rather thin and the inter¬ 
face may be disordered, a diffusive bulk treatment of the transport in the 
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normal metal is not generally applicable and the spin accumulation in the 
normal metal is in those situations not well defined. The appropriate vari¬ 
ables are then the torque on the ferromagnet, and the (electrical) current in 
the normal metal. The relations between current and torques are then found 
from symmetry considerations and microscopic calculations [SIIE21E3], and 
one finds that the general phenomenology remains the same. Moreover, the 
system is described phenomenologically in terms of an effective spin Hall 
angle and effective mixing conductance that characterize coupling between 
magnetic dynamics and current, and loss of angular momentum at the in¬ 
terface, respectively [33]. The drift-diffusion theory discussed here gives an 
expression for these parameters in terms of the bulk spin Hall angle 9sh, 
spin-mixing conductance and other parameters, that is appropriate if 

the diffusive treatment applies [33j . 

In summary, in this section we have established that temperature differ¬ 
ences and differences between magnon chemical potential and electron spin 
accumulation drive spin transport across interfaces between normal metals 
and magnetic insulators. Furthermore, we have demonstrated how a current 
through a conductor that is parallel to its interface with a magnetic insu¬ 
lator sets up a nonzero spin accumulation via the spin Hall effect, and how 
spin current injected from a magnetic insulator to normal metal through 
an interface can be detected via the inverse spin Hall effect in the normal 
metal. 

We emphasize that the above analysis shows that the metallic reservoirs 
effectively act as grand-canonical baths for the magnons in the magnetic in¬ 
sulator (see Fig. [It]), and that in this way the system we consider is rather 
close to the textbook grand-canonical treatment of Bose-Einstein condensa¬ 
tion. In the next section we discuss how by tuning the driving forces, i.e., 
electrical currents and temperature differences, a quasi-equilibrium Bose- 
Einstein condensate of magnons can be achieved and maintained. 


1.3 Pumping of quasi-equilibrium magnon condensation by 

spin current 

In the previous section we have reviewed recent developments in spintronics 
concerning the linearized interaction between electrons and thermal magnons 
at interfaces between magnetic insulators and normal metals. In this section 
we consider the more general situation of a magnetic insulator that is par¬ 
tially condensed, as a result of interactions with the metallic reservoirs. 

To make the discussion concrete, we consider a magnetic insulator that is 
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in the bulk described by a Heisenberg-model hamiltonian, 


^[S] = — 4 
^ ^ 2 h^ 


E s.-s,+5: 


<i,j> 


K 


2% 


q-Z 

2 


+ jSi. 

h 


( 1 . 8 ) 


with J the nearest-neighbor exchange energy, K > 0 the easy-plane anisotropy 
constant, and H > 0 the external field in units of energy. At zero tempera¬ 
ture and without dissipation, the dynamics of the average spin (Sj) ~ hSui 
is governed by the Landau-Lifshitz equation 


d{Si) 

dt 


1 

h 


{S^) X 


dH[{S)] 

d{S,) 


(1.9) 


which describes precessional dynamics around the effective field. 

The presence of the Bose-Einstein condensate is signalled by a nonzero 
expectation value T = (6) of the long-wavelength annihilation operator. 
Employing the same linearized Holstein-Primakoff transformation as in the 
previous section, we have that the magnetization direction for the condensed 
phase is at zero temperature given by n ~ (i/^/sReT, — A/2/sIm'k, no/s — 
1) where T = with no the condensate spin density and (p the 

azimuthal angle of the magnetization with the x-axis. In the homogeneous 
situation, the Landau-Lifshitz equation then results in 


h-^=/ic, (1-10) 

with the condensate chemical potential iJ.c = B — KS -|- KSno/s. The easy- 
plane anisotropy leads to a mean-held self-interaction for the condensate. 
Within the Landau-Lifshitz description, the condensate density is time in¬ 
dependent. We now incorporate spin injection at the interface and magne¬ 
tization relaxation to derive a rate equation for the condensate density. 


1.3.1 Condensate rate equation 

For simplicity, we consider the situation as depicted in Fig|1.2[ and focus only 


on the right reservoir. The interface spin current, determined by Eq. (1.1), 
is then given by 

27rs dt 


Using Eq. (1.10), and substituting Be ^ Be — BR to account for the nonzero 


spin accumulation in the right reservoir, this result is rewritten to yield the 
condensate contribution to the spin current across the interface between 
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metal and insulator m, 




( 1 . 11 ) 


The condensate spin current is accompanied by the spin current carried by 
thermal magnons, which, using Eq. (1.2), is found to be 

T^, Tr) = ^ f deD{e) (e - ^lR) 

’ vrs J 


riB 


e - 
ksTm 


, k-R\ 


( 1 . 12 ) 


In the above expressions for both condensate and thermal magnon spin cur¬ 
rents, the first term corresponds to spin pumping while the second accounts 
for spin transfer. 

Relaxation processes are at small energies accurately described by the 
Landau-Lifshitz-Gilbert phenomenology [Ml- This implies that the equa¬ 


tion for the dynamics of the magnetization direction in Eq. (1.9) acquires 
a Gilbert damping term parametrized by a dimensionless Gilbert damping 
constant a, given by 


dn 

dt 


dn 


= —an X 


rel 

Eor the condensate, this results in a loss of condensed magnons according to 


dno _ 2anoHc 

dt ~ h ' ^ ’ 

This loss term is similar in form as the contribution from spin pumping, as 
spin pumping corresponds to loss of angular momentum from the magnetic 
insulator to the normal metal across the interface. 

When the magnons are partially condensed, the density of thermal magnons 
is fixed and any spin current entering the magnetic insulator through its in¬ 
terface with the right reservoir is eventually absorbed by the condensate 
jl5j . We now assume that this absorption is instantaneous. Physically, this 
corresponds to the regime where the scattering rate due to interactions be¬ 
tween the thermal cloud and condensate is fast compared to the rate of 
magnon absorption/excitation at the interface, and to the damping rate. 
We denote jx = —jl^xif'ci fJ-R,Tm,TR) which is (minus) the interface spin 
current of thermal magnons in the partially condensed phase. The density 
of condensed magnons is then determined by the equation 


dnp _ jx _ gu^o ifJ'c - tJ-R.) _ 2 «rao/ic ^ jx _ no 

dt hLm 2 TTsLm tl fl ^Lm Tc ’ 


(1.14) 
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having defined 


1 

Tr 


1 

To 


g^lKSnp _ 2aKSno 

2TTs‘^hLm hs 


where tq denotes Tc in the limit np —>■ 0. For temperatures ksT K the 
thermal spin current jx is to a good approximation independent of no M- 
The four possible choices of absolute and relative signs of tq and jx give rise 
to four distinct quadrants in the steady-state {dnp/dt = 0) phase diagram 


(see Fig. 1.4): In region I both condensate and thermal spin current lead to 
loss of magnons and prevent formation of the condensate. In region II the 
condensate spin current leads to decay of condensate magnons, which are 
replenished from the thermal cloud by injection of thermal magnons. This 
results in a steady state in which a condensate exists. Region II crosses over 
to region III which is the swasing regime [3S], defined by the spin injected 
directly into the condensate as the spin accumulation is then larger than 
the magnon ground state energy. In this crossover the number of condensed 
magnons becomes larger, more rapidly so at lower temperatures. Finally, 
in region IV 2 the condensate spin current is not large enough to overcome 
the losses due to spin current by thermal magnons. In region IVi, on the 
other hand, two steady-state solutions exist. Only one of these is stable, 
however, and hence the number of condensate magnons needs to be above a 
critical number in order for the condensate to be maintained. This leads to 
hysteretic behavior in the first-order transitions from regions IV 2 and III to 
region IV 1 . The phase diagram in terms of the spin accumulation fiji and the 
magnon and electron temperatures is determined by calculating jx and Tp in 
terms of them, and reported in Ref. [2]. We conclude that quasi-equilibrium 
magnon Bose-Einstein condensation can be achieved by appropriate tuning 
of these driving forces. 


1.3.2 Spin superfluidity 

One of the most striking consequences of Bose-Einstein condensation of 
magnons is the resulting spin superfluidity |36[ [37] . Developing a theory for 
spin transport across the pumped magnon system is not straightforward, as 
the phase diagram itself will change in an inhomogeneous situation with re¬ 
spect to the homogeneous case discussed so far. Here, we ignore such issues 
and discuss spin superfluidity for the case of an equilibrium condensate. For 


the model hamiltonian in Eq. (1.8), an equilibrium condensate forms at low 


temperatures when the external field is small enough so that the easy-plane 
anisotropy tilts the magnetization away from the z-axis. The general phe- 
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Figure 1.4 Phase diagram of a quasi-equilibrium Bose-Einstein condensate 
maintained by spin-current injection. The point O corresponds to thermal 
equilibrium, i.e., no driving forces. The possible phases that meet at the 
critical point P are separared by transitions (solid lines) or cross-overs 
(dashed lines). For further details see main text. 


nomenology of the superfluid spin transport is expected to be similar when 
comparing equilibrium and pumped condensates. 

The spin superfluid transport is for a small precession angle and at low 
temperatures conveniently described by the continuity equation for the den¬ 
sity of condensed magnons (that includes a loss term due to Gilbert damp¬ 
ing), and the Josephson equation for the condensate phase. These are found 
from the Landau-Lifschitz-Gilbert equation [Eq. (1.9) with the Gilbert damp¬ 
ing added] and in the long wavelength limit given by (at zero temperature) 


dno 

h—= 
dt 

dvs _ 
dt 


-V • - 2anoiic 

2JgSI 
¥~ ’ 


(1.15) 


with Vg = —2JsVip/h the superfluid velocity and = riQhvs the condensate 
(superfluid) spin current. In the above we assume the condensate to be 
spatially homogeneous, and that no <C s and a <C 1. The stiffness Jg = 
JSa^j2, with a the lattice constant. For the configuration in Fig. 


1.2 


we 


find, by applying the boundary condition in Eq. (1.11) at both reservoirs. 
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that the spin current injected into the right reservoir is given by 


Js = 





g^l + 2TTasL 

m 


where we took hr = 0. Consequently, spin superfluidity is signalled by an 
algebraic decay of the spin current, and the ensuing inverse spin Hall voltage, 
as a function of the size of the magnetic insulator in the spin current direction 
[371 [Ml [39]. This is in contrast to the exponential decay (with the length 
scale im) that is found in the normal state m- The algebraic decay is a 
consequence of the small but nonzero Gilbert damping. In fact, the above 
shows that in the case when damping is absent (a = 0) the spin current 
through a spin superfluid magnetic insulator is limited only by the interface 
(via the mixing conductance 

For the geometry of our example, spin superfluidity may be pinned by 
anisotropies that break the rotation symmetry around the z-axis mm 
(which ultimately gives to the C/(l) symmetry that is spontaneously broken 
by the condensation). As an example, we co nsid er the addition of a term 
to the hamiltonian in Eq. (1.8). With this addition, the 


spin hamiltonian for constant density reduces to the energy 


E = dx. 


JsUo iV9f + {B- KS)no + 


KSnl 

2s 


— KxSno cos^ 6 


, (1.16) 


where we again took the long wavelength limit. The spin-current-carrying 
state requires gradients in the phase. These are penalized by the term ~ Kx 
which pins the phase at 0 = 0 or vr. The competition between exchange and 
in-plane anisotropy in the above energy thus defines a lower critical current 

jc,low — 271 o \/ JsSKx /, 


below which spin superfluid flow is pinned. Physically, this lower critical 
current follows from comparing the energy of the current-carrying state to 
the energy of a domain wall in 6 from, e.g., 0 to tt. Once the energy of the 
current-carrying state exceeds this domain wall energy, domain walls, and 
therefore gradients in the phase, are created, thus allowing finite superfluid 
spin currents. Below this lower critical current the anisotropy makes it en¬ 
ergetically favorable for the phase to remain homogeneous away from the 
interface [MISI]. 

The upper critical current is found by realizing that the condensate density 
is actually no = s{l — nl)/2 in terms of the z-component of the magneti¬ 
zation direction. From the above energy we then find that the superfluid 
flow is unstable towards increasing for |V0| ~ i? — KS, which defines 
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an upper critical current. When the current approaches this critical value, 
the superfluid-carrying state is relaxed by spontaneous vortex-induced phase 
slips transverse to our quasi-one-dimensional geometry m- 


1.4 Perspectives 

In this Chapter we have discussed how spin currents flowing across the inter¬ 
face between a magnetic insulator and normal metal can be used to achieve 
quasi-equilibrium Bose-Einstein condensation of magnons in the magnetic 
insulator, and, moreover, to probe its spin transport properties. Future theo¬ 
retical works should improve on our treatment of interactions, in particular 
regarding the interactions between condensate and thermal magnons and 
the role of phonons. Moreover, the description of the condensed phase may 
require inclusion of the dipolar interactions. Another direction for study are 
inhomogeneous situations and in particular the coupled spin-heat transport 
properties of the partially-condensed magnon system. 

On the experimental side, the best-studied systems are devices with only 
one metallic reservoir and where the magnetic insulator is YIG and the metal 
is Pt. More complicated set-ups involving more than one reservoir have not 
yet been seriously addressed. Moreover, from a materials-science perspective 
there is ample room for exploring novel materials and optimizing interface 
spin transport properties, e.g., by considering antiferromagnets [l2]. In con¬ 
clusion, we expect that the interplay between magnonic many-body physics 
and spintronics will be a source of new physics in the years to come. 
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the European Research Council (ERC) (R.A.D and A.B.), by the EU-FET 
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